Neutron Moderation in the Oklo Natural Reactor and the Time Variation of a 
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In previous analyses of the Oklo (Gabon) natural reactor to test for a possible time variation of 
the fine-structure constant a, a Maxwell-Boltzmann low energy neutron spectrum was assumed. We 
present here an analysis where a more realistic spectrum is employed and show that the most recent 
isotopic analysis of samples implies a decrease in a, over the last two billion years since the reactor 
was operating, of {apast — ce„ow)/ce > 4.5 x 10~* {6a confidence). Issues regarding the interpretation 
of the shifts of the low energy neutron absorption resonances are discussed. 

PACS numbers: 06.20.Jr, 28.41.-i, 28.20.Gd 



I. INTRODUCTION 

It was recognized by Shlyakhter Q that isotopic ra- 
tios of fission products and secondary neutron absorp- 
tion reactions that occurred with the Oklo natural reac- 
tor phenomenon some 2 billion years ago could be used 
to test whether the fine-structure constant a varies with 
time. Specifically, if a neutron absorber (e.g., ^''^Sm) 
has a low-energy neutron absorption resonance, and the 
excited compound state can decay by gamma emission, 
the resonance energy will change in a predictable way if 
a varies. The overlap of the resonance with the reac- 
tor neutron spectrum will subsequently change thereby 
altering the ratios of nearby isotopes. 

In addition to Q, two other detailed analyses of the 
Oklo phenomenon have been completed 0, 1^ ■ The pur- 
pose of this note is to assess the validity of the assumption 
in these analyses that the neutron spectrum can be de- 
scribed by a Maxwell-Boltzmann (MB) distribution, and 
to determine the effects on the limit of the time variation 
of a for different neutron spectra. The possibility of a 
non-MB spectrum was discussed in |3] (footnote 3), but 
was not considered to be important. However, for the 
high-accuracy results presented in Q, the effects due to 
deviations from a MB spectrum become very important. 



II. THE OKLO PHENOMENON AS A 
HOMOGENEOUS REACTOR 

The fact that the Oklo phenomenon occurred was in- 
dicated when uranium taken from this mine was found 
to be deficient in ^•^^U; subsequent tests confirmed this 
discovery. From core samples and knowledge of reactor 
dynamics, the age of the Oklo phenomenon was deter- 
mined to be about 2 billion years. The conditions neces- 
sary for the multiple reactor zones to exist can be used 
to place constraints upon many parameters such as con- 
centrations of various elements, temperature, moderator 
type and concentration, etc. while the reactor was oper- 
ating. For example, it can be shown that only hydrogen 



or deuterium could have served as moderators at Oklo 
and hydrogen present in water or as hydration of the ura- 
nium oxide appears to be the only plausible choice. As a 
first approximation to the reactor geometry, it seems rea- 
sonable to assume that the Oklo phenomenon was a ho- 
mogeneous reactor of infinite extent. This simplifies the 
estimation of the multiplication factor k and the neutron 
fiux ^{E, T) which we will see depends on the amount of 
hydrogen and neutron absorbing impurities in the ore. 



With the present isotopic ratio of ^^^U/^'^*U, the max- 
imum k that can be obtained in a homogeneous mixture 
of water and natural uranium is fc = .838, which oc- 
curs where there are 2.43 molecules of water per atom 
of uranium Q (see table 18.1, p. 612). The fractional 
isotopic abundance of ^ssy natural U at present is 
0.711%. Extrapolating this result implies a ^ssy relative 
isotopic concentration of about 1.2% to achieve k = 1. 
Two billion years ago, the relative concentration of ^ssy 
was about 3.7% (due to the different Hfetimes of ^ssy and 
^^^U). This implies that two billion years ago. A: « 1.3 for 
a system with 2.4 H2O molecules/U atom. If such a large 
fractional concentration of water was present, the Oklo 
reactor would have been highly divergent. Of the possi- 
ble effects that can reduce k, the two most likely are a low 
water (or hydrogen) concentration in the U deposit, and 
the presence of impurities. It is interesting to note that 
the reactor would have been somewhat self-stabilizing as 
a runaway to high temperatures would drive water out 
of the ore deposit. 



Precise knowledge of the nature and concentration of 
neutron-absorbing impurities at the time the Oklo re- 
actor was running is lacking. However, we can make a 
reasonable estimate of the possible maximum absorption 
by impurities. The forgoing arguments can made rigor- 
ous by estimating A; as a function of the hydrogen atomic 
number density fraction, — ph/pu and the atomic 
fraction of absorbing impurities fi = Pil pu, relative to 
the number density of U atoms. With ^35^ isotopic frac- 
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tion X, for thermal neutrons: 



pu [xaf 5 + (1 - x)gI^^] + pH(jH + P^(Jl 
2.47e--'^03i(i//ff)=** 

580x 



690x + 2.7(1 ~x) + + Y.^ /^< 



(1) 



where v = 2.47 is the number of neutrons released per 
^^^U fission, p is the resonance escape probability ob- 
tained from Eqs. (10.29) and (13.3) of 4], af^ = 580 b 
is the 235 u fission cross section, a'^^ — 690 b is the total 
2^^U total neutron absorption cross section, a'^^ = 2.7 
b is the 238u ^q^^^I neutron absorption cross section, 
(j^ = 0.33 b is the hydrogen absorption cross section. 
The total impurity absorption can be parameterized by 



/3 = E/»' 



(2) 



Furthermore, the effects of oxygen have been neglected. 
For X = .037, fc as a function of fn, and as a function 
of impurity concentration parameterized by /9, is shown 
in Fig. 1. With no impurities. A: = 1 is achieved when 
Jh ~ 1-5. For a higher hydrogen concentration, the reac- 
tor would have run away and the resulting high tempera- 
tures would have driven water out of the ore; the system 
was likely self-stabilizing to k — 1. As the impurity con- 
centration, parameterized by [3 in Fig. 1, increases, fu 
required for k ~ 1 increases also; a.t /3 = 5, fn — 2.6 
to achieve k = 1. For /3 > 11 b, it is not possible to 
achieve fc > 1. It is also known that the reactor oper- 
ated until X ~ 0.02 which implies a more restrictive limit 
of /3 < 3 b, for which fn = 5.8. We note that critical- 
ity calculations at these high 235, 238^ concentrations are 
subject to a number of corrections (see PP- 463-468) 
so the figures derived here are to be considered as pro- 
viding guidelines for the operation of the natural reactor 
which indicates that Jh is about 3, with (3 = 2. 

A ubiquitous substance that offers a significant thermal 
neutron absorption cross section is sodium chloride. The 
thermal absorption cross section for sodium is aa = 0.5 b, 
but for chlorine Ca = 33 b which limits fci < -1. The 
relative atomic concentrations of other impurities can be 
similarly estimated. 



III. CALCULATION OF & FOR A NON-MB 
NEUTRON SPECTRUM 

It is well-known that when neutrons are moderated 
in the presence of a 1/w absorber (such as U), the low- 
energy part of the spectrum becomes suppressed relative 
to a high-energy 1/E component of the flux. The relevant 
parameter that describes this suppression is (0, pp. 335- 
340) 



A 



2ASa(vT) 



(3) 



where 'Eaivx) = pa-o'aivr) is net absorption cross section 
of the homogeneous mixture for a neutron rms velocity 
corresponding to the moderator temperature, Eg = Ps-(Js 
is the net scattering cross section of the moderating 
atoms, and A is the atomic mass of the moderating 
atoms. In the present case, for neutron energies about 1 
eV, the atomic mass A = 1 for H which is responsible for 
most of the neutron moderation. When A sa 2, the 1/E 
tail contains the bulk of the neutrons. This is illustrated 
in Fig. 2 where the flux as a function of A/2 as obtained 
from a Monte Carlo calculation following the procedure 
in 5\. 

For the parameters derived in the last Section, that 
is Jh = 3, X = 0.037, and /3 = 2, A « 1. However, 
for neutron energies less than 1 eV, molecular effects be- 
come important; it is no longer a good approximation to 
assume that the neutron moderation occurs by scatter- 
ing from free atoms. The chemical bonding effect can be 
estimated by setting ^ = 2, so in the low neutron energy 
region (.001 < E < 0.1 eV), A « 2 at 300 K, and scales 
with temperature as 



(4) 



The elements ^'*^Sm and ^^^Gd have low energy neutron 
resonances and are particularly interesting to test for 
possible time variation of a. The formalism developed 
in Q, 1^ for the effective cross section of a substance 
with a low-energy neutron resonance can be understood 
as follows. The interaction of a neutron with a nucleus 
can be parameterized through the optical model, where 
the Hamiltonian is 



H = UiE) + iW{E) 



(5) 



where the absorption rate for neutrons with energy E is 
T{E) = 2W{E)/h. The absorption cross section can be 
expressed as 



a{E) 



V r 



fc2 iE-Er)yh^+Tt,//A 



(6) 



where hk is the neutron momentum, 5 is a parameter that 
depends upon the spins of the neutron and interacting 
nuclei, Er is the resonance energy and r„_^ are the rates 
for neutron and gamma ray absorption respectively. The 
parameters g, Er.and Tn,j for elements of interest here 
can be found in [3, 0] ■ 

The absorption cross section at a particular energy can 
also be described as the product of the rate of absorp- 
tion and the times t = L/v{E) spent in the material, 
normalized to the propagation length L, 



(t{E) = T{E)/v{E). 



(7) 



Far away from E^., r(i5) is nearly constant and (t{E) = 
aoVo/v{E), where Vq is the velocity at which the cross 
section was measured. Nuclei whose cross sections can 
be described this way are know as "1/u" absorbers. This 
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is a more efficient representation for thermal neutron ab- 
sorption of many nuclei such as U. 

The net absorption cross section is given by the integral 
of the cross section time the neutron flux ^{E), 



(a) = J a{E)^E)dE. 



(8) 



In assessing the effects of a varying low-energy resonance, 
because all 1/v cross sections follow the same universal 
function which have exactly the same form up to a con- 
stant, it is useful to define an effective relative cross sec- 
tion for the nuclei with a varying T{E) as 



/ a{E)^{E,T)dE 
J vop{E,T)dE 



(9) 



where cf{E) is the energy-dependent cross section, T is 
the moderator temperature, ^{E, T) and p{E, T) are the 
neutron flux and energy distributions, respectively, and 
vq is the velocity at which \/v absorption cross sections 
are measured. 

In this analysis and the others that have been per- 
formed so far, the finite temperature and motion of the 
absorbing nuclei is ignored. This is a small correction 
to the resonance width because A > 100 for the nuclei 
being considered, so the absorber velocity is much lower 
than the neutron velocity. Assuming A = 2 at 300 K, 
as a function of T, A = 2-\/300/T. A procedure similar 
to that described in W\ (numerical integration) was used 
to calculate the change in ij as a function of change in 
resonance energy, but replacing the MB spectrum with 
spectra derived from Fig. 2. The same resonance param- 
eters as employed in Q were used, and the numerical 
integration was tested with a MB spectrum for which 
the same results as Fig. 1 of were obtained. 

It is possible to determine the value of the cross section 
for ^^^Sm, (Ti49, from detailed isotopic analyses of core 
samples from Oklo and similarly for ^^^Gd to determine 
(T155. Such an analysis was performed in 3] and indicates 
that 



0-149 = 91 ± 6 kb. 



(10) 



This result can be used subsequently to determine the an- 
cient value of Er which can be compared to the present- 
day measured value. 

Results of the numerical integration for ^''^Sm are 
shown in Fig. 2 as a function of moderator tempera- 
ture with A = 2 at 300 K. In the mean temperature 
of the Oklo reactor was estimated to be about 600 K. 
From the curve for 600 K in Fig. 3 (expanded around 
the 91 ± 6 kb in the lower plot), we find that the change 
in resonance energy is 



AEr 



-45 



+7 
-15 



X 10-^ eV 



(11) 



which indicates more than a 6cr deviation from zero. Note 
that the "two solution" problem does not really exist in 
this case, but the upper magnitude range is expanded. 



In the calculation of the moderated neutron spectrum, 
effects of chemical bonds (which suppress the low energy 
spectrum) have been approximately included. As the 
atoms become more tightly bound, — > 00 so (A — > 00), 
and in this limit the neutron flux spectrum becomes pure 
l/E; the result of a calculation for this case is shown in 
Fig. 4. In this case, for which two solutions exist, 

AEr = (-135 ± 5) X 10"^ eV or (-58 ± 5) x 10"^ cV. 

(12) 

This result is roughly consistent with the A = 2 result; 
reality likely is somewhere between the A =2 and 00 
results. However, the result for A = 2 clearly indicates a 
non-zero change in Er at a level of more than 6a. 

Results for a calculation of ^^^Gd are shown in Fig. 5. 
The interpretation of the ^^'^Gd data is complicated by 
the presence of impurities, and an extensive analysis is 
given in 0. The Oklo data imply (7155 ~ 28 ± 10 kb, 
which, from Fig. 5, indicates, for the left-hand solution. 



AEr = (-90 ± 20) x lO"-' eV 



(13) 



This is consistent with the -'^^^Sm result because impuri- 
ties shift Er slightly, and the simple 1 /E spectrum over- 
estimates Er- 



IV. INTERPRETATION OF THE RESULTS 

A shift in resonance energy can be related to a shift in 
a by 



Aa AEr 

a 



Eo 



(14) 



A reasonable estimate of £'0 = — 1 MeV was presented in 
0, and is roughly the same for ^^^Sm and ^^^Gd. Taking 
the result Eq. (|11|) determines 



45 



-15 
-7 



X 10" 



(15) 



which is non-zero by more than 6(T. The central value 
plus twice the upper uncertainty can be used to set a 
95% confidence limit on the magnitude of the integrated 
rate of change over the last two billion years 



1-1 < 3.8 X 10-1 Vyr (95%conf) 



(16) 



By assuming higher temperatures, this limit decreases 
slightly. However, a significant modification of this re- 
sult requires temperatures beyond physically reasonable 
values as assumed in Q. 

It is expected that the neutron absorption resonance 
parameters (e.g., the neutron and gamma widths), other 
than Er, do not change significantly with a. In estimat- 
ing Ef) it was assumed that Er represents a fixed energy 
above the ground states of nuclei with A and A+1, and 
Eq is approximately the difference in the Coulomb en- 
ergies of the ground states. This result is in agreement 
with a more sophisticated analysis 
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V. CONCLUSION 

If a 1/v absorber is present in the moderator of a 
neutron chain reactor, the neutron energy spectrum can 
deviate significantly from the a Maxwell-Boltzmann dis- 
tribution. The Oklo phenomenon was rich in uranium 
which is a strong \/v absorber. By considering reaUstic 
deviations from a MaxweU-Boltzmann low-energy neu- 
tron spectrum, we have shown that a recent analysis of 
the Oklo natural reactor phenomenon implies an non- 
zero (> 6ct) variation of a over the last 2 billion years. 
The result is of opposite sign and an order of magnitude 
smaller than a recent astrophysical determination of the 
change in a 0, which is non-zero, but over a different 



time scale. It has been suggested that the time varia- 
tion of a is not monotonic ■ Our results are probably 
accurate to within 20%; a full MCNP model of the re- 
actor, assuming reasonable estimates of the temperature 
and impurities are possible, would be useful in light of 
the importance and interest in the possibility of a vary- 
ing a. Finally, these results might be interpreted more 
efficiently in terms of niq/ Kqcd for which the sensitivity 
is about two orders of magnitude higher than the sensi- 
tivity to a variation 0. 

We thank Tuan Nguyen for critical comments and in 
particular for clarifying the issues with the sign of the 
effect. 
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FIG. 2: Monte Carlo generated neutron velocity spectra for various values of A/2 = T,a{vT) /^s from 1 to .1, where vt ~ 
^/2kT/mn is the rms neutron velocity, assuming all moderation occurs from scattering on hydrogen, compared to a MB 
spectrum. 
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FIG. 3: Ca lculatio n of (T149 =< (J149 > as a function of change in resonance energy Er and as a function of temperature, 
for A f« 2v'300/r. Blue dots, 300 K; blue dot-dash: 400 K; blue dashes: 500 K; solid blue: 600 K; red dashes: 700 K; red 
dot-dash: 800 K. 




FIG. 5: ai5B for a MB spectrum (T =573 K) and a 1/E spectrum. 



